Mon. Not. R. Astron. Soc. 000, [Tj{5] (2009) Printed 18 December 2009 



(MN lATtfi style file v2.2) 



Spitzer IRAC imaging photometric study of the massive 
star forming region AFGL 437 

Lokesh Kumar Dewangan 1 *, & B.G. Anandarao 1 ! 

1 Astronomy & Astrophysics Division, Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India. 



ABSTRACT 

We present Spitzer IRAC mid-infrared photometry on the massive star forming region 
AFGL 437 (IRAS 03035+5819). From the IRAC colour-colour diagram, we identify 
several new embedded YSOs within 64 arcsec of the central compact cluster. Using the 
IRAC ratio images, we investigate the molecular outflows associated with the highly 
embedded young stellar object WK34 in the central cluster. We attribute the lobes seen 
(extended to ~ 0.16 pc in the north) in the ratio map to shocked molecular hydrogen 
emission. IRAC images reveal a large diffuse nebulosity associated with the central 
cluster and extending up to ~ 8.0 pc from south-west to north-east direction with its 
brightness gradually increasing from 3.6 to 8.0 /*m. A dense box-car-shaped nebula 
(more than 2.0 pc in size) situated to the south-west of the cluster shows molecular 
hydrogen emission that may have been caused by shock waves from the compact 
cluster sources. It seems that these sources are also responsible for the infrared-bright 
nebulosity. Using a 2D radiative transfer model, we derive from the spectral energy 
distributions, the mass, age and luminosity of all the YSOs identified within the central 
cluster. The SED modelling shows that the driving engine of the outflows, WK34, 
appears to be massive but very young and deeply embedded. The weighted mean 
values of the masses and ages of the 21 YSOs derived from the model are in the range 
1-10 M and io 41 ~ 6 - 4 yr respectively; while their luminosities are in the range of 

10 0.47-3.48 L 

Key words: stars: formation - stars: pre-main-sequence - infrared: ISM - ISM: jets 
and outflows - stars: winds and outflows 



1 INTRODUCTION 

Compact clusters embedded in giant molecular clouds pro- 
vide an opportunity to study recent star for mation over a 
wide r ange of masses in a small volume (e.g., lLada fc Ladal 
(2003)). Due to the large visual extinction suffered by the 
protostars in such clusters, infrared, especially mid- and far- 
infrared, observations help in studying them. AFGL 437 
(IRAS 03035+5819; G139.909+0.19 7) is one such compact 
embedded cluster of size ~ 15 arcsec (IKleinmann et al . 1977; 
lLada fe Ladal 120031) sit uated at a distance of 2.0 ± 0.5 kpc 
jArauilla fc Goldsmithll 19841) and having a tot al luminosity 



of AFGL 437N. The bright sou rces in the cluster su ffer a 
large extinction of A„ ~ 7 mag (|Cohen fc Kuhilll977l ). Ra- 
dio observations revealed that AFGL 437W and S, identi- 
fied as early B type ZAMS sta rs, are associated with Ultra- 



Compact (UC) H II r egions JWv nn- Williams et al 



Torrefies et"all Il992l: iKurtz, Churchwell fc Wood 
Wvnn- Williams et al.l l|l98ll ) and iTorrelles et al 



1981; 



1994) 



(|1992| ) 



of ~2 x 10 4 L Q (jWeintraub fc Kastnerl 1 199rj ). The cluster 



is associated with an optical reflection nebula and contains 
at least four embedded sou r ces called AFGL 437N, S, E 
and W JCohen fc Kuhill 19771). iRainer fc McLean! i| 19871 ) and 
later IWeintraub fc Kastnerl (1 19961 ) resolved a highly embed- 
ded source known as WK 34, ~2.5 arcsec to the South-East 
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have detected water masers towards AFGL 437N and W, 
whic h indicate th e yout h of the region of star forma- 
tion. iGomez et alj (Il992l ) found a compact, poorly col- 
limated bipolar CO outflow oriented in North-South di- 
rection. Studies of the clu ster in near-IR polarimetry 
(IWeintraub fc Kastnerl Il996h an d diff raction-limited imag- 
ing in 3.8 (jm (IWeintraub et al.lll996o have attributed the 
outflow to the highly-embedded, low-luminosity source WK 
34; this wa s later confirmed by the HST polarimetric imag- 
ing study dMeakin et al.1 120051 ). HST results showed fur- 
ther tha t the reflec t ion n ebulosity is due to the source 
WK 34. iDavis et al.l (1 19981 ) found FU emission 'wisps' to- 
wards the south-west of the compact cluster, but not associ- 
ated with the outflow. Recent high-resolution C 18 study 
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of the region by ISaito et all i|2006l . I2007T ) revealed a few 
dense cores/clumps around the central cluster in AFGL 
437, indicating the activity of star formation in the region. 
Ide Wit et all (|2009h provided spatially resolved 24.5 fim ob- 
servation on the main sources of the clu ster. Sub-mm/mm 
observations of the cluster were made bv lDent et aD (|l998l ) 
using beam sizes of 16-19 arcsec wh ich do not resolve the 
cluster members. iDevine et al, (2008) determined the age of 
the cluster to be 1-5 Myr. In the background of a number of 
these important observations and inferences made already 
on this interesting massive star forming region, our motiva- 
tion for the present study has been to look into the available 
(space) infrared imaging observations beyond 3 /xm, which 
are not studied so far, in order to understand the evolu- 
tionary stages of the highly embedded cluster members as 
well as the regions beyond the cluster and compare with the 
earlier results. 

Spitzer Infra-red Array Camera (IRAC) provides an op- 
portunity with an unprecedented high spatial resolution in 
thermal infrared wavelength regime that is very useful in 
identifying embedded sources in massive star forming re- 
gions. IRAC has four wavelength bands (with A e ///AA, 
3.55/0.75, 4.49/1.0, 5.73/1.43 and 7.87/2.91 H which in- 
clude molecular emissions such as those from H2 and Poly- 
cyclic Aromatic Hydrocarbon (PAH) molecules. The aims of 
the present study are to identify the embedded sources using 
the IRAC bands in order to classify the different stages of 
their evolution; to use ratio maps of the four IRAC bands 
in order to identify possible H2 em ission regions (o r PAH 
regi ons) following the s uggestion of ISmith fc Rosen! (|2005l ) 
and iPovich et all l|2007i ); and to present the nebulosity as- 
sociated with the cluster in all the IRAC bands. Further, 
we construct spectral energy distributions (SEDs) of the 
identified embedded sources in the cluster using not only 
the IRAC channels but also observations in J HK bands 
ilWeintraub fc Kastnerl Il996l ; iMeakin et aT]|2005h and mid - 



IR (|de Wit et al.ll2009l ) and sub-mm/mm ^Dent et alii 19981 ) 
regions; and model the S EDs using the on-line tool devel- 
oped bv lRobitaille et all (|2006l ). 

In Section 2, we describe the data used for the present 
study and the analysis tasks utilised. Section 3 presents the 
results on Spitzer IRAC photometry of 21 embedded sources 
associated with the dense young cluster. In the same section, 
we discuss our results in the light of what is already known 
about AFGL 437. Further in Section 3, we present the results 
and discussion on the ratio maps and on the SED modelling 
of the YSOs. In Section 4, we give the conclusions. 



2 SPITZER IRAC DATA ON AFGL 437 AND 
DATA REDUCTION 

The Spitzer Space Telescope (SST) IRAC archival images 
were obtained fr om the Spitzer pu blic archive, using 'leop- 
ard' software (see lFazio et al.l (|2004l ). for details on the IRAC 
instrument). The observations relevant for AFGL 437 were 
taken in High Dynamic Range (HDR) mode with 12s inte- 
gration time in all filters. These observations were a part 
of the project entitled, "The Role of Photo-dissociation Re- 
gions in High Mass Star Formation" (Program id 201; PI: 
G. Fazio). Basic Calibrated Data (BCDs) images were pro- 
cessed for 'jailbar' removal, saturation and 'muxbleed' cor- 



rection before making the final mosaic using Mopex and 
IDL softwares. A pixel ratio (defined as the ratio of the area 
formed by the original pixel scale, 1.22 arcsec/px, to that of 
the mosaiced pixel scale) of 2 was adopted for making the 
mosaic (giving a mosaic pixel scale of 0.86 arcsec/pixel) Q. 
Using these procedures, a total number of 60 BCD images 
of 5.2 x 5.2 arcmin 2 were mosaiced to make a final image of 
17.4 x 14.2 arcmin 2 in each of the four bands. Aperture pho- 
tometry was performed on the mosaic with 2.8 pixel aperture 
and sky annuli of 2.8 and 8.5 pixels using APPHOT task 
in IRAF package. The zero points for these apertures (in- 
cluding aperture corrections) were, 17.08, 17.30, 16.70 and 
15.88 mag for the 3.6, 4.5, 5.8, 8.0 /xm bands, here onwards 
called as Chi, Ch2, Ch3 and Ch4 respectively. The photo- 
metric uncertainties vary between 0.01 to 0.22 for the four 
channels, with Ch3 and Ch4 on the higher side. Ratio maps 
were produced from the IRAC channel images by using the 
standard procedure, with a pixel ratio of 8 that reveals the 
features prominently. 

For the purpose of SED modelling, the JHK photo - 
met ric data were t aken from IWeintraub fc Kastnerl (Il996l ) 
and IMeakin et al. I (I2005D as well as from 2MASS archives 
( Skrutskie et al ]|2006l). The 2 4.5fj,m photometric fluxes were 
obtained from Ide Wit et al.l |2009) (at a diffraction-limited 
resolution of 0.6 arcsec), for the three main sources, WK34, 
S and W, in the comp act central c luster . In addition, the 
sub-mm/mm data from lDent et alj (Il998l ) were used as up- 
per limits on the central sources, being unresolved in the 
large beam sizes. The details of the data compiled on in- 
dividual sources are given in section 3.3. We modelled the 
S EDs so constructed usin g the on-line SED-fitting tool due 
of iRobitaille et all l|2006l ). A criterion \ 2 ~ Xtest < 3 was 
chosen to obtain weighted means and standard deviations 
of individual physical parameters from sets of models for 
each object. For the SED modelling, we fo l lowed a similar 
procedure as described by [Grave fc Kumarl (|2009h . 



3 RESULTS AND DISCUSSION 

Fig 1 shows the IRAC images of the entire extent of AFGL 
437 in all four bands. The central compact cluster is shown 
marked by a square box in the Ch3 image (bottom left). 
As shown in Fig 1, the IRAC 8.0 /im band image (bottom 
right) reveals a diffuse bubble-/fan-like nebulous structure 
associated with and extended from the compact cluster in 
the south-west to the north-east direction (with a size of ~ 
8.0 pc). We find that brightness of the nebulosity gradually 
increases from 3.6 to 8.0 /im (Chl-Ch4), as seen in Fig 1. One 
can also notice in Fig 1, several filamentary structures prob- 
ably due to the UCHII regions associated with AFGL 437S 
and W in the compact cluster. Just below the box in the 
south-west in Fig 1 (between the two arrows in Ch3 image 
in bottom left), one can notice a dense cross-/boxcar-shaped 
structure that seems to be expanding into the surrounding 
interstellar medium (ISM). It has bright components in NE 
& SW (separated by ~ 1.70 pc) and in SE & NW (sepa- 
rated by ~ 1.5 pc) directions. The corners of this boxcar-like 
structure are seen faintly (as 'wisps') in the narrow band H2 
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image bv lDavis et al.l (| 19981 ). We find that this feature be- 
comes more prominent in longer wavelengths (see Fig 1). It 
may be possible that this was generated through interaction 
of stellar wind from the sources AFGL 437W & S with the 
surrounding dense material. The compact cluster (inside the 
box in Fig 1: Ch3) is shown enlarged in Fig 2 (in a colour 
composite of three images: Chi (blue), Ch2(green) and Ch4 
(red)), in which we have marked the central main sources 
namely AFGL 437 N, S, W and WK34 (E is not detected 
in IRAC) as well as other YSOs identified in the region (see 
Section 3.1). 

The compact cluster appears very nebulous in all 
the four bands. As s hown by the earlier authors (e.g., 
IWvnn- Williams et al.l l|l98ll )). the source AFGL 437W is 
associated with a blister HII region. It appears like a diffuse 
source with its brightness increasing progressively in the four 
IRAC channels Chl - 4 (an d is also seen in the 24.5 /im im- 
age bv lde Wit et al.l l)2009h ). A dense filamentary structure, 
seen to the right of AFGL 437W in all the bands, is prob- 
ably associated with the source and extends to about 38.3 
arcsec (0.4 pc) in NE-SW direction (see Section 3.2). 

3.1 Cluster Sources 

IRAC [3.6]-[4.5] vs [5.8]-[8.0] colour-colour diagram is shown 
to be a very powerful tool to classify proto-stellar objects 
into their evol utionary stages, s uch as Class 0/1, Class II 
and Class III (| Allen et al.ll2004l ). Such a colour-colour di- 
agram for the sources is shown in Fig [3] along with divi- 
sions for various pre-main-sequence classes, shown by boxes 
for Class 0/1, Class II and Class III sources. The criteria 
for Class I/II, which s hares the properties of Class I and 
Class II, are taken from lMegeath et al.l (|2004l ). Based on the 
above criteria, we identify 13-14 sources (including AFGL 
437S and WK34) as Class 0/1 and 4-5 as Class I/II includ- 
ing AFGL 437W, in the vicinity (within 90 arcsec) of the 
central compact cluster (see Fig [2}. Tabled] gives the photo- 
metric magnitudes of these sources, marked as si, s2 etc. in 
Fig [2] It may be noted that only 18 sources were detected 
in all the four channels. We have included in Table 1 three 
sources (s2, s3 and AFGL 437N) that were not detected in 
all the four bands, but satisfied one of the two colour crite- 
ria for Class I/II or Class II. In the Object column in Ta- 
ble 1, the numbers i n parentheses correspond to t he sources 
identified earlier bv lWeintraub fc Kastnerl l|l996h (see their 
Table 1 and F ig 2) . Some of these sources a re seen in the K 
band image of IWeintraub fc Kastnerl fl996), but are not de- 
tected in J and H bands. From IRAC photometry, we have 
identified a number of new YSOs, apart from the sources 
AFGL 43 7W, S, WK34, s2, s5-s8 and s lO that are classified 
earlier by We intraub fc Kastnerl |l996). A large number of 
these sources are highly embedded and suffer large extinc- 
tion. We have determined the visual extinction A„ for the 
sources that have JHK p hotometry l)Weintraub fc Kastnerl 
ll996l : lMeakin et al.l[2005h and found that WK34 is the most 
embedded of all, with A„ = 30-35 mag. In addition to the 
Class I and Class I/II cluster sources, about 35 Class II 
sources (from Fig 3) are also identified outside the cluster, 
spread over the nebular bubble. These are shown as open cir- 
cles in Chi image in Fig 1 (top left). Most of these sources 
occur within a distance of about 1.6 pc (165.4 arcsec) south- 
east of the cluster. For about 23 of these sources, good qual- 



ity 2MASS JHK photometric colours provide confirmation 
of their low mass Class II nature. This supports our con- 
jecture that the filamentary structures and the expanding 
boxcar-like nebula are possibly driven by the central cluster 



3.2 IRAC Ratio Maps 

As mentioned earlier, the IRAC bands contain a num- 
ber of prominent molecular lines/features. Chi contains H2 
vibrational-rotational lines while Ch2-4 mostly contain pure 
rotational lines. Chi, 3 and 4 also contain the PAH fea- 
tures at 3.3, 6.2, 7.7 an d 8.6 ym; but Ch2 do es not in- 
clude any PAH features. Ide Muizon et alj £l990) have de- 
tected some of the above mentioned PAH fea tures in AFGL 
437. As mentioned earlier, iDavis et al.l (|l998l ) have reported 
some 'wisp'-like features of H2 emission (in 2.121 /im narrow 
band filter that contains the 1-0S(1) line) at a few locations 
in AFGL 437. Several authors have utilised the ratios of 
IRAC bands to identi fy some of the molecu la r diagnostics 
menti o ned above (e.g., Smith fc Rosenl (|2005l ); |povich et"afl 
(|2007h ; lNeufeld fc Yuan I ( |2008T )L Since it is difficult to as- 
sess the contribution of different molecular transitions to 
different channels, the ratio maps are only indicative; un- 
til/unless supplemented by spectroscopic evidence. 

The Ch2 is more sensitive to H2 lines of high excita- 
tion temperatures while the Ch4 represents rotational lin es 
of low excitation temperatures (|Neufeld fc Yuan! {2008)). 
Likewise, the Ch2 does not have any PAH features while Ch4 
has. Thus, in the ratio image of Ch2/Ch4, the brighter re- 
gions indicate emission regions from higher excitations from 
H2 and the darker regions indicate PAH emission. This trend 
is reversed in the image of Ch4/Ch2 (i.e., bright regions show 
PAH emission and dark regions the H2 emission). 

Fig 4 shows the ratio maps of Ch2/Ch4 (a) and 
Ch4/Ch2 (b) in a region surrounding the central cluster. The 
bright regions in Fig 4a probably correspond to H2 emission; 
and those in Fig 4b correspond to PAH emission. The ratio 
map Ch2/Ch4 (Fig 4a) brings out well the cross-like fea- 
tures towards the south-west edge of the clust er; positions 
of whi ch match well with the faint 'wisps' seen in lDavis et al.l 
( 1998). These emission features could be the result of the ex- 
pansion of HII regions or stellar winds interacting with the 
local ISM. We have looked for PAH features by examining 
the ratio images of Ch4/Ch2; Ch3/Ch2 and Chl/Ch2 (since 
Ch2 does not have any PAH features). Fig 4b shows the ra- 
tio map of Ch4/Ch2. The ratio map shows bright regions 
around AFGL 437W, s5, s7 and slO; as well as towards north 
of the NE corner of the boxcar nebula. The ratio Ch3/Ch2 
(possible indicator for the 6.2 /im PAH feature) does not 
show any bright features. The requisite UV photon flux for 
exciting the PAH features possibly comes from the source 
AFGL 437W. The ratio map shows the boxcar/cross shaped 
feature in dark, probably because of lack of PAH emission. 
The bright narrow linear filament seen to the right of AFGL 
437W (between the dashed lines in Fig 4b) indicates prob- 
ably the ionization front fro m the massive star in the wake 
of which the PAH is excited dPovich et all (120071 ) ). The de- 
tection of the 'PAH filament' corroborates the blister mode l 
for AFGL 437W presented by I Wvnn- Williams et all l|l98ll ). 

We have examined more closely the ratio map of 
Ch2/Ch4 to look for the infrared counterpart of molecular 
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Table 1. Spitzer IRAC 4-channel photometry (in mag) of the YSOs identified in the central cluster of AFGL 437; the numbers in 
parentheses in the Object column are Weintraub & Kastner (1996) designations (see text) 



Object 


RA [2000] 


Dec [2000] 


Chi 


Err-Chl 


Ch2 


Err-Ch2 


Ch3 


Err-Ch3 


Ch4 


Err-Ch4 


class 


si (22) 


3:07:20.81 


58:30:35.60 


12.65 


0.03 


11.97 


0.01 


11.56 


0.19 


10.02 


0.12 


I 


s2 (23) 


3:07:22.53 


58:30:45.58 


9.65 


0.01 


8.97 


0.01 


7.97 


0.04 


— 


— 


— 


s3 (29) 


3:07:23.04 


58:30:48.90 


8.86 


0.04 


8.16 


0.03 


— 


— 


— 


— 


— 


W (28) 


3:07:23.88 


58:30:50.25 


8.34 


0.05 


7.77 


0.05 


5.50 


0.05 


3.52 


0.05 


1 


N (35) 


3:07:24.29 


58:30:54.84 


— 


— 


— 


— 


5.23 


0.04 


3.34 


0.04 


— 


S (20) 


3:07:24.53 


58:30:42.90 


6.72 


0.01 


5.46 


0.01 


4.11 


0.00 


2.68 


0.01 


1 


WK34 (34) 


3:07:24.55 


58:30:52.76 


7.24 


0.03 


6.05 


0.01 


4.98 


0.03 


3.51 


0.05 


1 


s4 


3:07:25.94 


58:30:08.90 


12.36 


0.02 


10.55 


0.01 


9.26 


0.01 


8.35 


0.04 


I 


s5 (7) 


3:07:26.21 


58:30:20.96 


9.71 


0.02 


9.54 


0.02 


6.52 


0.02 


4.73 


0.02 


I/II 


s6 (31) 


3:07:26.45 


58:30:52.62 


9.47 


0.01 


8.89 


0.02 


6.63 


0.01 


4.53 


0.01 


I 


s7 (11) 


3:07:26.50 


58:30:25.42 


9.60 


0.02 


9.40 


0.02 


6.63 


0.02 


4.85 


0.05 


I/II 


s8 (36) 


3:07:26.53 


58:31:08.32 


10.76 


0.03 


10.24 


0.02 


8.98 


0.02 


7.32 


0.04 


1 


s9 


3:07:27.38 


58:30:12.28 


11.97 


0.02 


11.16 


0.01 


10.51 


0.06 


9.36 


0.09 


I 


slO (40) 


3:07:27.40 


58:31:15.76 


10.05 


0.02 


9.66 


0.02 


8.36 


0.06 


6.46 


0.06 


I/II 


sll (17) 


3:07:27.76 


58:30:35.32 


10.66 


0.02 


9.69 


0.01 


8.87 


0.04 


7.66 


0.05 


1 


sl2 


3:07:28.69 


58:30:47.39 


12.68 


0.03 


11.22 


0.01 


10.81 


0.08 


10.30 


0.22 


1 


sl3 


3:07:30.22 


58:30:58.88 


13.43 


0.02 


12.65 


0.02 


12.06 


0.09 


10.81 


0.12 


I 


sl4 


3:07:30.83 


58:31:37.84 


15.62 


0.07 


14.03 


0.02 


12.99 


0.14 


11.49 


0.16 


1 


sl5 


3:07:31.32 


58:31:12.45 


13.20 


0.01 


12.53 


0.01 


11.74 


0.03 


10.35 


0.03 


I 


sl6 


3:07:31.56 


58:29:59.03 


11.86 


0.01 


11.67 


0.01 


9.22 


0.01 


7.18 


0.01 


I/II 


sl7 


3:07:32.60 


58:31:23.03 


13.94 


0.02 


13.72 


0.02 


12.70 


0.13 


11.11 


0.22 


I/II 



outflow originating from the highly embedded YSO WK34. 
Fig 5 shows a zoomed-in image of Ch2/Ch4 ratio contours 
over-laid on the Ch4 image. One can notice the prominent 
lobe/outflow stretched northwards of WK34 and slightly 
bent towards NE, with a total extension of ~ 0.16 pc. The 
outflow direction and size are consistent with the earlier re- 
ports. The bending itself is attributed to the presence of t he 
nearby source AFGL 437N (|Weintraub fc Kastnerl (|l996h ). 



3.3 SED modelling of Cluster Sources 

By modelling the SEDs of the identified YSOs, we derive 
physical parameters of both their photospheres and mass 
accreting disks/envelopes. For this purpose, we constructed 
the SEDs for all the 21 sources listed in Table 1, from the 
optical to sub-mm data depending upon the availability in 
the archives and published literature. The optical datqj are 
available only for the sources W (B,V bands), s7 (V upper 
limit) and sl6 (B,V,R upper limits). In fact, the lack of opti- 
cal data is indicative of the large extinction that most of the 
sources suffer. The JHK fluxes for the sources slO, sl6 and 
s!7 are taken from 2MASS a nd for s2, s3, s5-s8 and sll from 
I Weintraub fc Kastnerl ill99rj). While the source si has only 
K band flux from Weintraub fc Kastnerl l|l996j ). s4, s9, and 
sl2-sl5 do not have near-infrared counterpa rts. The JHK 
photo metric fluxes for WK 34 are taken from iMeakin et all 
1211115 :. We have used sub-mm/mm data (|Dent et all (|l998l )) 
as upper limits for WK34, S, W and N (because of the 
large beam sizes that nearly encompass the entire central 
region of the cluster), only to constrain the models. The 
rest of the sources listed in Table 1 are outside the field 
of view of the sub-mm observations. The SEDs are then 
modelled using an on-line 2D-radiative transfer tool due of 



from http:/ /vizier. u-strasbg.fr/viz-bin/VizieR 



iRobitaille et al.l (|2006l ) , which assumes an accretion scenario 
with a central star surrounded by an accretion disk, an in- 
falling flattened envel ope and bipolar cavitie s. This tool was 
successfully tested bv lRobitaille et ail (|2007|). on a sample of 
low mass YSOs and bv lGrave fc Kumar (2009) on high mass 
protostars. The SED model tool requires a minimum of three 
data points that are of good quality. In our sample all the 
21 sources have at least three such data points. For a source 
that has meager number of data points, clearly the tool picks 
out a large number of solutions that can fit the data well, 
within the specified limit on \ 2 - R the SED has larger data, 
spread over the wavelength region of 0.5 to 1000 /im, then 
the model would be better constrained to yield r esults with 
least standard deviations l|Robitaille et al.| [2007). The dis- 
tance to the source and visual extinction are to be given as 
input parameters, usually as a range of values. This leads to 
a further degeneracy of the models. In our case, however, the 
cluster distance is fairly well-determined and the extinctions 
were estimated from the available JHK data on individual 
sources. But, in order to avoid 'over-interpretation' of SEDs 
l|Robitaille et al.ll2007h . we have provided a range of visual 
extinction values for each object to account for the uncer- 
tainties in their determination from JHK photometry. Only 
those models are accepted which satisfy the criterion of x 2 

- Xbest < 3. 

We show in Fig [6] the SED modelling results for the 
sources AFGL 437W, S, N and WK34. Table [2] gives the 
weighted mean values of the derived physical parameters for 
the four sources along with the standard deviations. The ta- 
ble also lists the model-derived weighted mean values of A„ 
with standard deviations and the degeneracy of the mod- 
els (i.e., the number of solutions that satisfy the criterion 
mentioned above). The model-derived parameters listed in 
Table 2 indicate that all the four sources are likely to be 
massive (early B type). This is consistent (within the stan- 
dard deviations) with the observations on W and S that are 
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Table 2. Physical Parameters derived from SED modelling of the main sources AFGL 437 W,S,N and WK34 (see text) 



Source 


Age 


M, 


T, 






Mdisk 




Degeneracy/ 


Name 


log(yr) 


M 


log(T(K)) 


log(L Q ) 


log(M yr- 1 ) 


log(M yr" 1 ) 


mag 


No. of models 


W 


5.36±0.19 


8.28±1.48 


4.28±0.11 


3.55±0.23 


_4.84±0.17 


-7.01±1.16 


4.0±0.7 


16 


S 


4.35±0.90 


7.79±1.72 


3.86±0.31 


3.14±0.39 


-4.12±1.30 


-5.84±1.12 


9.2±1.8 


72 


N 


4.82±0.46 


9.36±1.62 


3.99±0.19 


3.48±0.30 


-3.96±0.94 


-6.02±0.84 


18.0±8.0 


132 


WK34 


3.94±0.57 


7.23±1.91 


3.71±0.10 


3.00±0.30 


-4.06±0.55 


-4.69±0.39 


23.0±5.6 


3 



associated with UCHII regions. It may be noted that obser- 
vationally very little is known about the spectral type of the 
source N. 

The model parameters for WK34 suggest that the 
source is massive but of young age with effective temper- 
ature still not sufficient to create an HII region (see Table 
2). Its low luminosity is suggested by earlier workers also. 
This is also reflected in its outflow, which shows emission 
in H2 (as revealed by the ratio map in Fig 5) but not in 
PAH; indicating that WK34 is still not hot enough to pro- 
duce sufficient UV flux. The modelling suggests that the 
source AFGL 437W and S are also massive stars but more 
evolved to attain su fficiently large ef fective temperature to 
excite an HII region. iQin et all l|2008h derived a CO outflow 
entrainment rate of 7.4xl0 -4 M /yr for the outflow source 
(WK 34), which agre e s with that observed for massive stars 
(e.g., iBeuther et ail (|2002h l. In compari son, the l o w ma ss 
stars show much less entrainment rates l|Wu et a l. (1996)). 
In a study of th e mole cular outflows from high mass YSOs, 
iRidge fe Moord (|200ll ) have concluded that these outflows 
are often poorly collimated, which again points to the pos- 
sibility that WK34 is a massive protostellar object. 

From the SED modelling of the 21 sources, we infer 
that the weighted mean values of the masses and ages of the 
YSOs in the central cluster (for the 21 sources listed in Table 
1) are in the ranges 1-10 M and io 41-6 - 4 yr respectively; 
while the luminosities are in the range of iq - 47-3 - 48 L . 



4 CONCLUSIONS 

The important conclusions of this work are as follows: 

(i) Spitzer IRAC imaging photometry is presented on the 
massive star forming region AFGL 437; 

(ii) Several new embedded YSOs are identified within 90 
arcsec of the central compact cluster; 

(iii) The IRAC ratio maps indicate molecular outflow cor- 
responding to WK34 which is possibly due to H2; 

(iv) The SED modelling of the outflow driving source 
WK34 indicates that it is a massive but very young pro- 
tostar not yet able to drive a HII region; 

(v) SED modelling of the 21 sources gives their masses 
in the range 1-10 M , ages io 41-6 - 4 yr and luminosities 

10 0.47-3.48 t 
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Right Ascension (J2000.0) 

Figure 1. IRAC images of AFGL 437 (size ~ 17.4 X 14.2 arcmin 2 ) are shown in all four channels (in log scale). The central compact 
cluster is marked by the square box (of side ~ 129 arcsec or 1.2 pc) in Ch3 image (bottom left). The associated bubble-like extended 
nebulosity is of size ~ 8.0 pc in the SW-NE direction. The cross-like or boxcar-shaped dense nebula expanding into the local ISM is 
shown marked by the two arrows in the Ch3 image (bottom left). The open circles in Chi (top left) mark the Class II sources identified 
from IRAC photometry. 
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03 h 07 m 30 s 25 s 20 s 

Right Ascension (J2000.0) 

Figure 2. The compact cluster (see the box in Fig 1) is shown in a colour composite zoomed-in image (8.5 (red), 4.5 (green) and 
(blue) (im). The identified YSOs are shown as si, s2.. along with the brighter sources AFGL 437S, W and WK34 (see Table 1). 



8 L.K. Dewangan & B.G. Anandarao 



in 
i 

to 1 



2.5 



2.0 



1.5 



1.0 



0.5 



0.0 



-0.5 



-1.0 



Class 0/I 



A 
A 

A «WK34 




A 

ClassJI ^ 



oo 



o <> 



pt> o 2 



w. . 



• t • 

* Class l/ll 



Class III 



_l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L_ 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

[5.8]-[8.0] 



Figure 3. Mid-IR colour-colour diagra m constructed usin g the Spitzer IRAC bands. The boxes indicate possible regions of Class 0/1, 
Class II, and Class III sources based on l| Allen et aLl fetKM). Class I/II sources share the properties of both Class I and Class II sources. 
The filled circles represent the sources within 64 arcsec of the centre of the cluster (see Table 1 and Fig 2); the asterisks are Class I/II 
stars and the diamonds are for Class II (see Fig 1 for their locations). The Class III stars are shown as plus-signs. The triangles are 
Class I sources found quite far away from the central cluster. The principal sources of the cluster, AFGL 437 S, W, and WK34 can be 
seen marked inside the box designated for Cl ass 0/1 protostars. Ext inction vector for K band, An =5 is shown in the diagram, which is 
calculated using averaged reddening law from iFIahertv et a.1.1 d2007t) ■ 
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Figure 4. IRAC Ch2/Ch4 (a) and Ch4/Ch2 (b) ratio images in log scale (from mosaics made with pixel ratio of 2). The central cluster 
and the cross- or boxcar-like structure to its south-west can be noticed. This structure has a size of 1.7 pc in the NE-SW and 1.5 pc 
in NW-SE directions. The bright patches in the image on the left (a) show regions dominated by H2 emission; while in the image on 
the right (b) the bright patches indicate regions that are emitting PAH bands. The prominent outflow lobe, some of the YSOs, and the 
H2 emission patches in the four corners of the cross-like structure are marked (by dashed curves). The near-horizontal tracks across the 
image arc artifacts. The filament-like structure to the right of AFGL 437W is marked by dashed lines in the upper half of the figure. 



10 L.K. Dewangan & B.G. Anandarao 



AFGL 437 



o 
o 
o 

CNI 

—> 

c 

_o 

-I— » 

o 

o 

CD 
Q 



20" 



10" - 



31'0" - 



50" - 



40" - 



58°30'30" 




03 h 07 m 28 s 



26 s 24 s 
Right Ascension (J2000.0) 



22 s 



Figure 5. IRAC 8 £tm zoomed-in image of the cluster is shown in log scale overlaid by the contours of IRAC ratio image of (Ch2/Ch4). 
Contours show the outflow lobes associated with the source WK34. Outflow lobes are stretched in north direction with a bend towards 
east. The IRAC Ch 4 (8 (im) image is made with a pixel ratio of 8 in a spatial extent of ~ 55 X 55 arcscc 2 ; while the ratio contours are 
plotted with a minimum of 0.178 and maximum of 8 M Jy/Sr. 
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Figure 6. Spectral Energy Distributions for four YSOs in the central cluster: AFGL 437W, S (top panels) N and WK34 (bottom panels). 
Filled circles are observed fluxes of good quality (with filled triangles as upper limits) taken from archives or published literature (see 
text for references) and the curves show the model fits. The thin black curve corresponds to the best model. The dashed curves represent 
photosphcric contributions. The model parameters of the YSOs are listed in Table 2. 



